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ABSTRACT 



The proper scientific analysis of a large amount of stellar 
spectra requires certain capabilities of the analysing tool 
(e.g. precise semi-automatic fitting for normalisation of 
the continuum or line-list assisted measurement of spec- 
tral lines) as well as flexible list-driven datafile handling. 

While most astronomical legacy packages comprise pow- 
erful analysing and data management features allowing 
rapid processing of quite complex data, current VO tools 
are lacking support of even basic capabilities conmionly 
used in stellar spectroscopy. 

Especially the high resolution optical spectroscopy of 
stars with rapid line profile variations or of those with 
complicated emission profiles benefits from a number of 
specific methods unsupported in todays Virtual Observa- 
tory (VO) tools and even lacking definition in VO spectral 
protocols. 

In our contribution we identify these techniques, describe 
their possible implementations and finally give a short 
overview of several VO-compatible tools emphasising 
their deficiencies and comparing with the capabilities of 
common legacy packages. 

Key words: spectroscopy, stellar continuum, spectral 
lines, radial velocity, web services, SSAP , Virtual Ob- 
servatory. 



1. INTRODUCTION 



There is a wealth of calibrated astronomical spectra ac- 
cessible in current data archives, however most of them 
are not suitable for direct physical analysis. Being pre- 
reduced by some automatic pipeline or individually by 
manual reduction, they are mostly stored in archives in 
form of FITS file or ASCII table as the relation of in- 
tensity in arbitrary numbers (instrumental counts, digital 
numbers etc.) and wavelength (or frequency). 



The scientific analysis of such spectra requires further 
processing by the variety of different methods. In certain 
studies a huge number of spectra has to be collected from 
different servers (e.g. in different spectral regions ) and 
transformed into common units. Often the transformed 
spectra or variables measured in them (e.g. equivalent 
width) are visually inspected as the time series ordered 
by time of observation or by the orbital phase. 

The tools with built-in VO protocols may considerably 
increase the efficiency of astronomical analysis by elim- 
inating the data collection, aggregation and conversion 
and by allowing the astronomer to fully concentrate on 
the visual properties of data (e.g. searching the variability 
of line profiles or constructing the radial velocity curve). 
To be accepted quickly by the wide astronomical com- 
munity a number of methods and recipes has to be im- 
plemented in a manner very similar to the current legacy 
applications. We attempted to identify the most useful 
methods and give hints for their implementation in VO 
clients. 



2. COMMON METHODS OF ASTRONOMICAL 
SPECTROSCOPIC ANALYSIS 



2.1. Rectification 

Usually the first step in data analysis is the rectifica- 
tion of the spectrum (normalisation of the continuum to 
level 1.0). Although absolutely flux-calibrated spectrum 
is common for space instruments, it is rather exceptional 
in ground-based optical stellar spectroscopy (used mainly 
for determination of spectral type from low resolution 
spectra). 

Usually, simple polynomials of low degree are applied for 
this purpose but there is a number of problems that have 
to be addressed: 

Anchoring of edges - all fitting polynomials have to be 
correctly anchored at the edges to work well. How- 
ever in many cases the edge of the chip has a dif- 
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ferent behaviour than the middle. Usually the edges 
show rapid decline in intensity due to the vignetting 
of optics or due to lower detector quantum efficiency 
on low resolution spectra covering whole visible 
range both at UV and IR regions. Such curve with 
flat middle and rapidly falling edges is very difficult 
to fit with low order polynomials. 

Flexibility of fitting curve - The big problem of fitting 
curves is connected with their flexibility. The curve 
should be rigid enough to bridge consistently the 
continuum gap over strong wide lines (Hydrogen, 
Helium) and flexible enough the fit well the changes 
in curvature of continuum. This is especially impor- 
tant for normalising spectra of hot stars (strong He 
and Hydrogen B aimer lines) stars or rapidly rotat- 
ing stars with broadened lines. Too high orders of fit 
tend to enter into the line, low order do not copy the 
tilt of continua at the edges well. 

Majority of legacy spectral tools use Chebyschev or 
Legendre polynomials and cubic splines (e.g. IRAF 
splot, MIDAS XAlice or Starlink DIPSO). What 
is needed is a fitting routine that mimics the freehand 
drawing of the experienced astronomer. One of the 
little known but very powerful procedure is INTEP 
piill, 1982), which uses Hermittean polynomials. It 
has been successfully used many yea rs for normal 



ising spectra of Be stars in SPEFO ( Skoda 



1996) 



and recently implemented in SPLAT-VO. Special at- 
tention was given to freehand drawing algorithms 
by |Akim a (1970). His procedure is implemented in 
SPLAT-VO as well. 

Range of data points included into fit Some programs 
(e.g. IRAF splot) allow to fit only the existing data 
points but some regions (e.g. wide lines) may be ex- 
cluded by giving the number of distinct regions. The 
more flexible approach (e.g. in SPEFO, SPLAT-VO) 
is the possibility of placing individual control points 
everywhere (e.g. above the apparent continuum). 
The extrapolation of fitting line beyond the available 
data range can help to decide about the best parame- 
ters (polynomial order, number of spline segments) 
from its global behaviour (smoothness of curve and 
asymptotic approach to the particular value). 

The problem of physical continuum The continuum 
may be not seen on some objects so it cannot be 
fitted using the available data points. Problem in 
cooler stars is to locate the real position of the 
continuum - it may not be present at all due to 
severe blend of lines or molecular bands - then a 
model is needed. 

It is extremely difficult to normalise certain stars 
with unusual profiles — like stars with rapidly 
expanding envelopes (post-AGB outbursts, nova 
shells), having strong emission lines with P-Cyg 
profile dominating whole spectrum. 

For some stars with extremely strong emission lines 
(e.g. some Be and B[e] stars) the dynamic range of 
16-bit ADC is insufficient to obtain the unsaturated 



line together with a well exposed continuum. Often 
the continuum is hidden in background noise. 

Normalising echelle spectra A really challenging prob- 
lem (which had not yet been solved fully) is the nor- 
malisation of echelle spectra - especially with wide 
absorption lines or complicated profiles (P Cyg, Be 
stars). Some wide lines may span several echelle 
orders and so the individual orders have to be pre- 
cisely unblazed and then merged. This problem 
is very complex, depending on construction of in- 
divi dual spectrograph, the obser v ation strategy etc. 
See 'Skoda & Hensberge" (2003'), 'S koda & §le"ctita 
([2004 ), Prugniel & Soubiran (2001 J or |Erspamer & 
[North| ( [2QQ2l ): ' 

fitting procedures The more advanced way of nor- 
malising is the least squares fitting based on some 
theoretical model. It is very popular in radioastron- 
omy and X-ray astronomy ( XSPEC ) rather for es- 
timation of the global SED fit parameters (usually 
only as combination of several simple power law 
functions). Quite useful in optical spectroscopy is 
the usage of appropriate synthetic spectrum to find 
the real continuum in data. The implementation of 
such a procedure in VO tools could call the The- 
ory VO (TVO) servers to get correct models (e.g. 
Kurucz) and perhaps to improve the models by re- 
computing on the GRID in a iterative way. Unfor- 
tunately for many interesting objects models (even 
spectral classification) are only speculative. 



2.2. Simple visual comparison 

A lot of the information about the behaviour of astronom- 
ical objects can be estimated just by a visual inspection of 
a spectrum (spectral type, peculiarity, emission) or a time 
series of spectra (pulsations, binarity). The basic method 
is the overplotting of many spectra in the same units and 
scale. It may be very efficient with VO-enabled tools ob- 
taining the number of spectra cached immediately from 
VO spectral servers through Simple Spectra Access Pro- 
tocol (SSAP). There are many possibilities of plots: 



• Overplotting spectra of different objects in the same 
region 

• Overplotting the same object in different wavelength 
ranges (from X or UV over IR to radio). If the spec- 
tra are flux calibrated, it is the construction of spec- 
tra energy distribution (SED). 

• Comparison of observation with models (rotation- 
ally broadened theoretical spectra) 

• Plot of different lines overplotted in radial velocity 
scale may give he information about physical struc- 
ture (e.g. for study of kinetic behaviour of differ- 
ent elements in expanding envelopes ( [Kipper et al.[ 
[2004j)). 
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Figure 1. Dynamical spectrum of Hj3 line profile vari- 
ability of 59 Cyg. Residuals from average profile of 38 
spectra are binned to 20 phase bins corresponding to pe- 
riod 28.192 days. Two expanded cycles shown for clarity. 
Individual profiles are overplotted above. After Maintz 
( |20g?l ) 



2.3. Study of line profile variations 

The high resolution spectra with high SNR may reveal 
on some objects small variations of the profile of spec- 
tral lines. The study of LPV requires many (even hun- 
dreds) of spectra to be overplotted to see the changes 
well. Sometimes the animation of changes in individual 
spectral line is very impressive. The low SNR data should 
be removed not to spoil the image. To identify them some 
form of colour highlighting should be implemented to- 
gether with single- stroke delete command to remove it 
from plotting list. Such a plot is helpful in asteroseismol- 
ogy or for estimating changes in stellar winds. 



2.4. Dynamical spectrum 

It is sometimes called the gray representation or trailed 
spectrum. The basic idea is to find the small time- 
dependent deviations of individual line profiles from 
some average. 

First the average of many high dispersion high SNR spec- 
tra (with removal of outliers) is prepared (called tem- 
plate spectrum). Then each individual spectrum in time 
series is either divided by the template (quotient spec- 
trum) or the template is subtracted from it (the differen- 
tial spectrum). The group of similar resulting intensities 
is given the same colour or level of gray. See Fig. [T] and 
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Figure 2. Dynamical spectrum of Silll line profile vari- 
ability ofX Sco. The average profile of 8 observing nights 
was used as a template shown below the panels. After 
\Uytterhoeven\^^004) 
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Fig.[2| More examples m ay be found e.g. in|de Jong et al. 
Maintz] ( [2003] ) or |Uytterhoeven| ( |2004| ). 



The result is drawn in 2D image where on horizontal axis 
is a wavelength in the line profile or corresponding RV 
relative to laboratory wavelength, on the vertical axis the 
time of middle of observation (in HJD) or the circular 
phase when the data are folded with certain period. Valu- 
able are interactive features, like zoom of whole dataset, 
removal of bad spectra from series, adjustable colour cuts 
and look-up-tables etc. 



2.5. Measurement of radial velocity and higher mo- 
ments of line profile 

The one of the important information received from spec- 
trum is the radial velocity of the object - e.g. the bi- 
narity can be revealed , or the possession of extrasolar 
planet. RV is sometimes presented as zero-th moment 
of the spectral line profile. The higher moments of the 
line profile are important as well. The first moment is 
equivalent width. The combination of higher moments of 
line profile is a one of the possible ways of determination 
ofnon radial pulsation modes numbers /, m (Aerts et al.[ 
[1992]). 



RV by fitting profiles This is a most common method 
implemented in every spectral package. The data 
points in given range of line profile are directly fit- 
ted by Gauss, Lorentz or Voight function. At least 
rough normalisation is required. Problems occur in 
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Figure 3. The estimation of RV by mirroring of profile. 
This case shows the match of line wings. 



case of asymmetric and complicated profiles (e.g. P- 
Cyg type). 

RV estimation by profile mirroring It is a method used 
already many years ago in oscilloscopic compara- 
tors. The goal is to estimate the best match of the 
line profile and its flipped image (mirror around y 
axis in certain wavelength). The mirrored profile is 
being shifted interactively until the best match is as- 
sessed. This method (as implemented in SPEFO) 
has been successfully used for measurement of RV 
of stars with very complicated emission profiles (e.g. 
Be stars, luminous blue variables with strong wind, 
nova outbursts). 

The implementation should allow the region of in- 
terest (wing or core) in vertical and the width of 
flipped area in horizontal direction to be easily ad- 
justable. The reference line position is best taken 
from given linelist providing rest wavelength. The 
strength of this method lies in the capability of mea- 
suring the match of core and match of wings sep- 
arately. See Fig. [3) More examples are given in 



Parimucha & Skoda (2006). This allows to study a 
different physics of particular line-forming regions 
(e.g. shells, jets, stellar winds). 

There is a close connection with bisector analysis, 
when the region of match is controlled in different 
levels of depth of the line profile. 



RV by cross correlation It is often used to measure the 
RV of late type stars with known spectral classifi- 
cation. It does not measure RV of individual 9iu 
lines but a RV shift of the whole spectrum. Many 
sharp lines are required and a good template (usually 
a synthetic spectrum convolved to the similar res- 
olution and rotationa l broadening as observed one) 
( [Barton et al.| |2000| ). Continuum normalisation is 
critical. 

Many legacy applications (IRAF package xrv and 
RVSAO, standalone program daospec, MIDAS task 
XCORRELATE . . . ). Often used for searching extraso- 
lar planets (cross-correlation with Iodine cell spec- 
trum). The cross-correlation is an ideal method for 
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Figure 4. The construction of bisector of line profil e (left 
panel) and the zoomed bisector ( right panel). From \Gray\ 
( [7952] ). 



measurement of RV in echelle spectra (many lines 
in different spectral regions may be combined to- 
gether). 



2.6. Measurement of equivalent width 

The measurement of equivalent width (EW) of a given 
line needs the determination of area bordered by line pro- 
file and the continuum. This area is transformed to the 
rectangular line with the depth of 1.0. The width of this 
rectangle is called equivalent width. Emission lines by 
definition have negative EW. 

The measurement is sensitive to continuum placement. 
The extremely shallow lines give large error in EW 
( [Vrancken et al.|p^97] ). EW are often used in abundance 
analysis. A good check of correct normalisation is the 
comparison of EW of spectral lines of the same element 
in different spectral regions ( [Hensberge et al.[[2000| ). 



The comparison of EW measured in spectrum of point- 
like solar system bodies (minor planets) with high res- 
olution solar spectrum from solar telescopes gives the 
information about the intrinsic spectrograph properties 
(stray light) an d quality of pipeline reduction ( ^Erspamerj 
[& North! [2002] ). 



2.7. Bisector analysis 

It is a method describing quantitatively the tiny asymme- 
try or subtle changes in line profiles. It is easily done by 
marking the middle of horizontal cuts of profile in differ- 
ent line depths. The line connecting such points (called 
bisector) is then zoomed in horizontal (wavelength or 
RV) direction. See Fig.|4] 

The characteristic shape of bisector gives the information 
about turbulence fields (e.g. convection) in stellar photo- 
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sphere , characterised by t he value of micro-turbulent ve- 
locity ( [Gray! p^82l |2QQ5 ) or about other processes caus- 
ing the tiny profile asymmetry. It has been used suc- 



cessfully for searching of extrasolar planets (Poyichet 
[aLj [2001 ) or in asteroseismology. Requires high reso- 
lution (echelle) and high SNR normalised spectra. Still 
some smoothing of edges in steeply declining line wings 
is helpful ( jMartinez Fiorenzano et al. 2005J). 



2.8. Rotational broadening 

Before comparing the synthetic spectrum with observed, 
the rotational broadening has to be applied, especially for 
rapid rotators. The estimate of rotational velocity v s'mi 
is thus obtained. 

Although simple formula of convolution with parabolic 
kernel has been mostly used, the precise computation is 
very complicated (the problem of limb darkenin g) espe- 
cially in non-LTE case ( [Hadrava & Kubat[|2003| ). 

Sometimes a Fourier transformation of a line profile (with 
high SNR) is applied to get the first estimate of rotational 
velocity ( |Gavrilovic & Jankov,,2007j . 



well for analysis of multi-pe riodical non-radial pul- 
sations (|de Jong et al.[|1999 ). 
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Figure 5. Example of 2D periodo gram of 811114552 A l ine 
profile variability ofX Sco. After Uytterhoeven\ ((2004^ . 



3. COMPLEX PROCESSING METHODS 



2.9. Period Analysis 

One of the most popular methods used in astronomy is 
the period analysis. Its aim is to find the hidden peri- 
ods of variability of given object. Sometimes this period 
can be identified with some physical mechanism (e.g. or- 
bital period of binaries, rotational modulation or pulsa- 
tions). Wide range of objects show the multi-periodicity 
on various time scales (e.g. binary with pulsating com- 
ponents). Ones the suspected period is found, the data 
may be folded accordingly, plotted in circular phase cor- 
responding this period. Very helpful is a interactive ca- 
pability of showing data folded while selecting different 
peaks at the periodogram. 

The number of methods applied is large, but there are 
several very popular in astronomy: 



• Folding techniques: Phase dispersion minimisation 
PDM ( ,Stellingwerf| [19781 ). 



• Fourier methods: Widely used algorithm CLEAN 
( Roberts et al. 1987|) or its modified Lomb-Scargle 
version ( |Lombl]T?76[[Scargle|p^89^ 



• Periodogram of RV of bisector (sometimes called lo- 
cal RV). The periods in position of bisectors in dif- 
ferent line depths may give an idea where the pro- 
file is affecte d by given non radial pulsation mode 
(|Koubsky"eral.,2007). 



• Periodogram of line profile. It is a kind of gray 
representation with period on vertical axis and the 
zoomed line profile on horizontal. See Fig.|5] Works 



For the complex techniques given below, a lot of addi- 
tional information is required in addition to spectral data. 
Programs require complicated configuration files in given 
format and some interactive trials to find the best results 
using the output from recent run as input to next one. 
They are written often in FORTRAN without graphical 
interface or even the plotting capabilities. They are de- 
signed for batch runs driven by parameter files. Their 
application is limited to special cases of certain objects 
with appropriate physical properties. 



3.1. Doppler imaging 



It was introduced by |Vogt & Penrod ( 1983a| ) as a method 
allowing the surface mapping of stellar spots. First test 
were done on st ars of RS CVn type and on ( Oph ( Vogt &] 
Penrod[p^983b| ). Works well on rapid rotators and needs 
a high resolution spectra with very high SNR (300-500). 
Tho whole rotational period should be covered well, bet- 
ter several times. When all the requirements are met, the 
map of surface features (spots, nodes of non radial pulsa- 
tions) is obtained with very high accuracy. 



3.2. Doppler tomography 



It was introduced by |Marsh & Horne ( |1988 ) for mapping 
the distribution of emitting circumstellar matter in binary 
system. One of successful application gave a picture of 
accretion jets in Algols (Richards, 2004). It uses trailed 
spectrum in velocity scale. The result is 2D image in ve- 
locity space. The transformation of radial velocity space 
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to coordinate space is ambiguous, which causes problems 
in interpretation of Doppler tomograms. 



3.3. Spectra disentangling 

This method allows to disentangle the spectra of individ- 
ual stars in binary or multiple systems even in case of 
heavy blending of lines. It supposes the changes in line 
profile are caused only by combination of Doppler shifted 
components (no intrinsic variability of star). The best so- 
lution of orbital parameters and disentangled line profiles 
of individual stellar components are found by least square 
global minimisation. Automatically removes the telluric 
lines with great precision. Requires good orbital coverage 
and the estimate of orbital parameters. Two approaches 
exist: 



• Blackbody fit 

• No RV measurement built-in 

• No complex operations with spectra allowed 

• It can be called directly from VizieR 

• Can work with linelists through Simple Line Access 
Protocol (SLAP) 

• Theoretical VO (TVO) supported (synthetic spectra: 
Kurucz, disks) 

• Rapid development 

• Integration of some methods in users own program 
possible 

• JAVA applet 

• Contains PLASTIC VO-interoperability layer 

• Dereddening of extragalactic objects built-in 

• Support of dimensional equations for units 
(DIMEQ, SCALEQ) 

• Rather complicated view of users data - needs to cre- 
ate SSA wrapper and prepare VOTable even for one 
simple FITS spectrum before viewing 



Wavelength space disentangling developed by |Bagn- 



uolo & Gies| ( |1991| ) and improved by |Simon & 
Sturm (1994]). It needs a large memory to store 
sparse matrices, requires large computing power. 
Gives more reliable results 



Fourier space disentangling introduced by [Hadrava 
(1995 ) and Hadrava, (,1997 ) in program KOREL. 
Another program available today (still based on 
KOREL ideas) is FDBINARY plijic et al.| |20Q4| ). 
They work in Fourier space, and transform the 
wavelengths into In A. They solve a small amount 
of linear equations, so they are memory savvy and 
can be run on even small computer. The method, 
however, requires perfect continuum fit (difficult to 
achieve at merged echelle spectra). On unreliably 
rectified spectra strange artifacts may appear plijic 
[etar||2QQtl ). 



4. VO-ENABLED TOOLS 



Although there is a number of various tools written 
for the analysis of spectra of astronomical objects (e.g 
IRAF, MIDAS, Starlink packages splot, spectool, 
DIP SO, XALICE as well as numerous single-purpose C, 
FORTRAN and IDL routines), the number of modern 
tools understanding the VO protocols (especially spectra 
access protocol - SSAP) is still very low. We are giving a 
brief overview of their capabilities and status of develop- 
ment at the time of the submission of these proceedings 
(May 2007). 



4.1. VOSpec 



• Current version 2.5 



http://esavo.esa.int/vospec/ 
Developed at ESA 

Very simple (mainly for building SED) 
Polynomial and Gaussian fits only 



4.2. SpecView 



Current version 2.13 

http://www.stsci.edu/resources/software_hardware/specview] 
It IS available as JAVA applet or stand-alone applica- 
tion 

Developed and supported by STScI 
Understands a number of formats from HST instru- 
ments and most NASA/ESA satellites 
Understands general FITS in form of binary tables 
Does not handle simple ID FITS (binary images 
with CDELT, CRVAL) 

Easy work is with local data (if binary or ASCII ta- 
bles) 

Well tailored for practical spectral analysis 

Powerful linelists (many included) 

Fitting of models by — even user models may be 

constructed 

Built-in number of standard stars spectra and models 
for various physical conditions 
Whole Kurucz library of spectra available for imme- 
diate display. 

Only simple polynomials for visual fit available 
Dereddening, CLOUDY models built-in 
Proper error propagation of fitted values 
Support of dimensional equations 
Does not support PLASTIC 



4.3. SPLAT-VO 

• Current version 3.7 

• http://star- www.dur.ac.uk/~pdraper/splat/splat- vo/| 

• Supported by JCM'l alter closing Starlink 

• Most advanced for stellar astronomy (fitting by free- 
hand drawing, INTEP) 

• RV measurement - both Gaussian fit and mirroring 

• Custom line list for individual spectra available 

• Wavelet analysis 

• Full featured data and editor and spreadsheet 

• Publication quality output, powerful plotting op- 
tions, annotations 

• Supported PLASTIC 

• Reads ID FITS image files (e.g. rebinned IRAF 
multispec files) 

• Asynchronous SSAP queries 
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4.4. Period04 

• Current version 1.0.2 

• http://www.univie.ac.at/tops7Penod047l 

• Its predecessor was widely used Period98 

• Period04 is rewritten in Java (and partly C) 

• Needs formatted text files 

• VO-interface not supported, only local VOTable 
files 

• Several period finding methods built-in 

• Can handle even the regular period shifts 



4.5. FROG 

• Part of Starlink Java package 

• Included in last JCMT Starlink release (Hokulei) 

• http://www.jach.hawaii.e du/software/starlink/| 

• VO-protocols built-in 

• Almost same capabilities as Period04 

• Unclear status of development after closing Starlink 
- all the links to documentation and files are now 
wrong 

• Last version since 2004 - no further development no- 
ticed 

• The period searching engine available remotely as a 
SOAP service 

• Web service for Fourier transform 

• Collaborates with TOPCAT, but not PLASTIC 

• Easy but powerful for doing period analysis in VO 
environment 



5. IMPLEMENTATION STRATEGIES 



As was shown above, the current VO-enabled tools are 
already very flexible to allow the large part of analy- 
sis of astronomical spectra to be directly performed in 
the VO environment using the data through SSAP. How- 
ever, many "classical " astronomers, working mostly with 
middle and high resolution stellar spectra in visual wave- 
lengths, will not use these tools until most of the features 
described earlier are implemented. 

The functional behaviour of such implementations is very 
important as well. To be quickly adopted by the con- 
servative part of stellar community, the tool should be 
self-learning: it should allow interactive work with one 
or couple of spectra to adjust the parameters (ranges for 
measurement, degree of fitting polynomials etc.). But af- 
ter tuning its behaviour it should be able to apply the same 
recipes to large number of spectra with similar character- 
istics. The very promising way how to do it is the concept 
of workflows. 

There is a important question where to put the engine for 
spectral analysis. All tools today re the JAVA applications 
running as VO-client on astronomer's computer. Adding 
more features (new methods of data processing) requires 
the new version of client to be downloaded and installed. 
Despite the possible automatic installation through Java 



WebStart there is a serious problem of size and complex- 
ity of such tool. Although the new functions could be 
added as client plugins (similar to Firefox plugins), we 
do not consider such a solution to be prospective on a 
long-term scale as it is against the original concept of 
VO. VO has been usually advertised as a facility allowing 
seamless federation of all distributed data archives and 
computing resources providing data processing services 
through Web Services. 

Certain amount of processing power should be thus left at 
the data providing servers. It may perform some simple 
operations before sending the spectrum to the client. For 
example it may transform the units (both wavelength and 
flux), cut out only the requested wavelengths (e.g. 20 A 
long section centred around Ha) or even the server might 
perform the normalisation of spectrum by approximate 
continuum fit - either computed using given parameters 
of fit or just applying the fit included in calibrated spec- 
trum (some reduction programs produce spectra in format 
containing both the unrectified spectrum with intensity in 
ADU and the continuum fit described by type and degree 
of polynomials and list of control points). 

But most of complicated processing (e.g. spectra disen- 
tangling, Doppler tomography, period analysis, should be 
rather performed on dedicated problem-oriented servers 
using Web Services or similar technology. It is, how- 
ever, difficult to convince the authors of powerful legacy 
(mostly FORTRAN) applications to re-implement them 
according to requirements of VO infrastructure (they do 
not know the JAVA, VO protocols etc.) and it is dan- 
gerous to ask the software developer to blindly re-code 
such a applications without detailed knowledge of the 
employed mathematical methods (including their bottle- 
necks and week points). 

Thus the only feasible way how to create most of such 
Web Services is to create the VO-compatible wrappers 
calling the legacy application as it was designed using 
its native input and output file format. The collaboration 
with authors is important anyway. 

The utilisation of dedicated servers allows easy extension 
of processing capabilities by adding another service to the 
same server (e.g. by adding Apache-like server module) 
or by installing new dedicated server. The natural way 
of solving complicated computing demanding problems 
is to deploy such services on GRID. 

There are still tasks in astronomical spectra analysis re- 
quiring visual interaction (e.g. precise fitting of compli- 
cated continuum, measurement of RV by profile mirror- 
ing or adjusting the colour scale of dynamical spectrum to 
strengthen particular feature in it). Such (usually simple) 
procedures may be left on the client side. Today, when 
there is a very few applications written in form of Web 
Services (to be deployed on server side), the extension of 
client capabilities seems to be the easiest solution even 
for some server-oriented tasks (e.g. the cutout of selected 
spectral region). 



8 



6. CONCLUSIONS 



Astronomical spectroscopy uses a wide range of tech- 
niques with different level of complexity to achieve its fi- 
nal goal — to estimate the most precise and reliable infor- 
mation about celestial objects. The large part of spectro- 
scopic analysis today has been accomplished by several 
independent non VO-compatible legacy packages, where 
each works with different local files in its own data for- 
mat. Analysis of large number of spectra is thus very 
tedious work requiring good data bookkeeping. 

Accomplishing the analysis in VO infrastructure may 
benefit from automatic aggregation of distributed archive 
resources (e.g. the multispectral research), seamless on- 
the-fly data conversion, common interoperability of all 
tools (using PLASTIC protocol) and powerful graphical 
visualisation of measured and derived quantities (e.g. in 
VOPlot or Mirage). 

By introduction of modern VO-aware tools into the astro- 
nomical spectral analysis a remarkable increase of effec- 
tivity of astronomical research can be achieved. 
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